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1 Introduction 



In February 2001, the Muon (g-2) Collaboration of the E821 experiment at the Brookhaven AGS 
released a new value of the anomalous magnetic moment of the muon , measured with an unprece- 
dented accuracy of 1.3 ppm [parts per million]. This annoucement has caused quite some excitement 
in the particle physics community. Indeed, this experimental value was claimed to show a deviation 
of 2.6 a with one of the most accurate evaluations of the anomalous magnetic moment of the muon 
within the standard model. It was subsequently shown that a sign error in one of the theoretical 
contributions was responsible for a sizeable part of this discrepancy, which eventually only amounted 
to 1.6 a. However, this event had the merit to draw the attention to the fact that low energy but 
high precision experiments represent real potentialities, complementary to the high energy accelerator 
programs, for evidencing possible new degrees of freedom, supersymmetry or whatever else, beyond 
those described by the standard model of electromagnetic, weak, and strong interactions. 

Clearly, in order for theory to match such an accurate measurement, calculations in the standard model 
have to be pushed to their very limits. The difficulty is not only one of having to compute higher orders 
in perturbation theory, but also to correctly take into account strong interaction contributions involving 
low-energy scales, where non perturbative effects are important, and which therefore represent a real 
theoretical challenge. Furthermore, in the meantime the members of experimental team at Brookhaven 
have further improved upon their measurement. In July 2002, they have announced a new result, in 
perfect agreement with the one they had obtained about one year and a half before, but with an error 
lowered down to the 0.7 ppm level. In addition, the theoretical evaluation of a specific and important 
contribution, called hadronic vacuum polarization, on which we shall have more to say later on, has 
shown a discrepancy between different data sets that are used as imputs. Depending on the choice 
between these conflicting data, the discrepancy between the experimental and the theoretical values 
can be as small as 1 a, certainly not a case for beyond the standard model physics, or can reach the 
3 a level, a much more promising situation as far as the possibility of "new physics" is concerned. 

The purpose of this account is to give an overview of the main features of the theoretical calculations 
that have been done in order to obtain accurate predictions for the anomalous magnetic moment of 
the muon within the standard model. Actually, all three charged leptons, e ± , fi^, and r^, of the 
standard model can be treated on the same footing, except that the very different values of their 
masses will induce different sensitivities with respect to the mass scales involved in the higher order 
quantum corrections. Thus, the anomalous magnetic moment of the electron is almost only [but not 
quite] sensitive to the electromagnetic interactions of the leptons, and its value is barely affected by 
strong interaction effects or by weak interaction corrections. On the other hand, the strengths of the 
latter two types of corrections are enhanced by a considerable factor, ~ (m T /m e ) 2 ~ 1.2 x 10 7 , in the 
case of the r, as compared to the electron. The same huge enhancement factor would also affect the 
contributions coming from degrees of freedom beyond the standard model, so that the measurement of 
the anomalous magnetic moment of the r would represent the best opportunity to detect new physics. 
Unfortunately, the very short lifetime of the r lepton which, precisely because of its high mass, can 
also decay into hadronic states, makes such a measurement impossible at present. The muon lies 
somewhat in the intermediate range of mass scales 2 and its lifetime still makes a measurement of its 
anomalous magnetic moment possible. However, in order to obtain an accurate [that is, comparable to 
the experimental accuracy] prediction, the contributions of all the sectors of the standard model have 

2 The corresponding enhancement factor is ~ (m M /m e ) 2 ~4x 10 4 . 
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to be known very precisely. Therefore, although the other two lepton flavours will also be discussed, 
since this does not really require additional work, the emphasis of these lectures will nevertheless be 
put on the muon. 

There exist several excellent reviews and introductions, which the interested reader may consult. As 
far as the situation up to 1990 is concerned, the collection of articles published in Ref. [1] offers a 
wealth of information, on both theory and experiment. Very useful accounts of earlier theoretical work 
are presented in Refs. [2, 3]. Among the more recent reviews, Refs. [4, 5, 6, 7, 8] are most informative. 
I shall not touch the subject of the study of new physics scenarios which might offer an explanation 
for a possible deviation between the standard model prediction of the magnetic moment of the muon 
and its experimental value, should such a deviation be confirmed in the future. For this aspect, I refer 
the reader to [9] and to the articles quoted therein, or to [10] for a list of the recent papers on the 
subject. 



2 General considerations 

In the context of relativistic quantum mechanics, the interaction of a pointlike spin one-half particle of 
charge and mass mi [£ stands hereafter for any of the three charged lepton flavours e, fx or r] with 
an external electromagnetic field A^x) is described by the Dirac equation with the minimal coupling 
prescription, 

in ~^ = ca ■ (^-ihV - ^Aj + (3m E c 2 + e e A tp . (2.1) 

In the non relativistic limit, this reduces to the Pauli equation for the two-component spinor p de- 
scribing the large components of the Dirac spinor tp, 



ih 



dp 



(-ihV - {e f Jc)Af 
2rri£ 



2mec 



er B + eeAo 



p. 



(2.2) 



As is well known, this equation amounts to associate with the particle's spin a magnetic moment 



M s = g e 



(—) 

\2mecJ 



S, S = %■ 



(2.3) 



with a gyromagnetic ratio predicted to be gi = 2. 



In the context of quantum field theory, the response to an external electromagnetic field is described 
by the matrix element of the electromagnetic current 3 J' P [spin projections and Dirac indices of the 
spinors are not written explicitly] 

<r ( P ')\J p (o)\r (p)> = u(p')rV,p)u(p) , (2.4) 

with [kf, =p'^- Pp] 

r p (p', P ) = Fi(k 2 )Y + -1- F 2 (k 2 )a^k v - F 3 (k 2 ) l5 a^k u + F±{k 2 )[k 2 Y - 2m^F] 7 5 ■ (2.5) 

3 In the standard model, J p denotes the total electromagnetic current, with the contributions of all the charged 
elementary fields in presence, leptons, quarks, electroweak gauge bosons,... 
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This expression of the matrix element {£~ (p')\J p (0)\£~ (p)) is the most general that follows from 
Lorentz invariance, the Dirac equation for the two spinors, (p 1 — m)u(p) = 0, u(p')(p'' — m) = 0, and 
the conservation of the electromagnetic current, (d ■ J) (x) = 0. The two first form factors, F\(k 2 ) 
and F2(k 2 ), are known as the Dirac form factor and the Pauli form factor, respectively. Since the 
electric charge operator Q is given, in units of the charge ee, by 

Q = J dxj (x°,x), (2.6) 

the form factor F\ (k 2 ) is normalized by the condition F\ (0) = 1. The presence of the form factor F3 (k 2 ) 
requires both parity and time reversal invariance to be broken, whereas F^k 2 ) can be different from 
zero provided parity is broken. Both F^(k 2 ) and F^{k 2 ) are therefore absent if only electromagnetic 
and strong interactions are considered [we leave aside the possibility of having a non vanishing vacuum 
angle in the strong interaction sector]. On the other hand, in the standard model, the weak interactions 
violate both parity and time reversal symmetry, so that they actually induce non vanishing expressions 
for these form factors. 



The above form factors are defined for k 2 < 0, and they are real in this region if the current J p {x) is 
hermitian. Due to general properties of quantum field theory, like causality, analyticity, and crossing 
symmetry, these form factors can be analytically continued into the whole complex k 2 plane with a 
cut for k 2 > Am 2 . They then become complex functions, obeying the Schwartz reflection property 
Fi{k 2 )* = Fi{k 2 *). For k 2 > Am 2 , the form factors Fi{k 2 + ie) describe the crossed channel matrix 
element {£~(p')£ + (p)\J p (0)\0). Furthermore, at k 2 = 0, they describe the residue of the s-channel 
photon pole in the S-matrix element for elastic £ + £~ scattering. 



At tree level in the standard model, one finds 

Ff rcc {k 2 ) = 1 , F,f rcc {k 2 ) = 0, i = 2, 3, 4 . 



(2.7) 



In order to obtain non zero values for ^(/c 2 ), ^(A; 2 ), and F^k 2 ) already at tree level, the interaction of 
the Dirac field with the photon field would have to depart from the minimal coupling prescription. 
For instance, the modification [T^ v = d^A v — d u A p , J p = ipj p ip] 



J d 4 xC int = ~J d 4 xJ p A p 

-> J d 4 x£ iQt 



I d 4 x 



with 



leads to 



J P = J P - ^—a e d p (ipaf,pip) - —d^iipi^appif)) , 



F^ik 2 ) = 1 , F 2 trcc (^) = a e , F 3 tree (fc 2 ) = di/et , F 4 tree (fe 



i tree ( i2 



tree / 7,2 a 



i tree / ;2 \ 



0. 



(2.8) 



(2.9) 



(2.10) 



4 The current J p is still a conserved four-vector, therefore the matrix element (£ (p')\J' p (0)\£ (p)} also takes the 
form (2.4), (2.5), with appropriate form factors Fi(k 2 ). 
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The equation satisfied by the Dirac spinor ip then reads 



ca ■ ^—ih'V — —Aj + (3m^c 2 + eiAo 



2rri£ 



a e p (ia ■ E - £ • B) - M e (3 (S • E + ia. ■ B) 



(2.11) 



and the corresponding non relativistic limit becomes 



ih 



dip 
~dt 



(-ihV - (e e /c)A) 2 e £ h 
2rri£ 2m£C 



(1 + ai)cr • B - hdecr ■ E + e e A + ■■■ 



(2.12) 



Thus the coupling constant a,£ induces a shift in the gyromagnetic factor, gt = 2(1 + a^), while dt 
gives rise to an electric dipole moment. The modification (2.8) of the interaction with the photon field 
introduces two arbitrary constants, and both terms produce a non renormalizable interaction. Non 
constant values of the form factors could be generated at tree level upon introducing [11] additional non 
renormalizable couplings, involving derivatives of the external field of the type U n A p ^ which preserve 
the gauge invariance of the corresponding field equation satisfied by ip. In a similar way, one can also 
introduce terms which induce a nonzero value for F4. In a renormalizable framework, like QED or the 
standard model, calculable non vanishing values for F2(k 2 ), Fs(k 2 ), and F±{k 2 ) are generated by the 
loop corrections. In particular, the latter will likewise induce an anomalous magnetic moment 
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(g e - 2) = F 2 (0) 



(2.13) 



and an electric dipole moment de = e^i^O); -F^O), which corresponds to an axial radius of the lepton, 
is also called the anapole moment [12, 13, 14], and is sensitive to the gradients of the external fields. 



If we consider only the electromagnetic and the strong interactions, the current J p is gauge invariant, 
and the two form factors that remain in that case, F\{k 2 ) and ^(/c 2 ), do not depend on the gauges 
chosen in order to quantize the photon and the gluon gauge fields. This is no longer the case if the 
weak interactions are included as well, since J p now transforms in a non trivial way under a weak 
gauge transformation, and the corresponding form factors in general depend on the gauge choices. As 
we have already mentioned above, the zero momentum transfer values -F«(0), i = 1,2,3,4 describe a 
physical S-matrix element. To the extent that the perturbative S-matrix of the standard model does 
not depend on the gauge fixing parameters to any order of the renormalized perturbation expansion, 
the quantities -Fj(O) should define bona fide gauge-fixing independent observables. 

The computation of T p (p',p) is often a tedious task, especially if higher loop contributions are con- 
sidered. It is therefore useful to concentrate the efforts on computing the form factor of interest, e.g. 
F 2 (k 2 ) in the case of the anomalous magnetic moment. This can be achieved upon projecting out the 
different form factors [15, 16] using the following general expression 6 

Fi(k 2 ) = tr [A P (p',p)(p" + m e )T p (p',p)^+m £ )] , (2.14) 

5 Terms involving the gradients of the external fields E and B or terms nonlinear in these fields are not shown. 
6 From now on, I most of the time use the system of units where fi = 1, c = 1. Other projectors on F 2 (k 2 ) have also 
been devised, see e.g. [17], but are not currently used. 
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with 



ZK ' k 2 k 2 -Amj k 2 {k 2 - Am}) 2 

AS(P',P) = -T^ ,, ' 2 -BV- (2.15) 



For fc — ► 0, one has 



and 



1 r kf r>P U i 

{p' + m i )k p 2 {p',p){p" + m,) = - y + mt )\- - + W - 2-)lL + ...]. (2.17) 
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The last expression behaves as ~ 1/A; as the external photon four momentum fc^ vanishes, so that one 
may worry about the finiteness of i*2(0) obtained upon using Eq. (2.14). This problem is solved by 
the fact that T p (p',p) satisfies the Ward identity 

(p'-p) p TP(p',p) = 0, (2.18) 

following from the conservation of the electromagnetic current. Therefore, the identity 

rV,p) = -k<r-^r°(p', P ) (2.19) 

provides the additional power of k which ensures a finite result as k^ — > 0. 

The presence of three different interactions in the standard model naturally leads one to consider the 
following decomposition of the anomalous magnetic moment af. 

QED . had , weak /r> orA 

ag = af + a,£ + dp . (2.2V) 

The first term, a^ ED , denotes all the contributions which arise from loops involving only virtual 
photons and leptons. Among these, it is useful to distinguish those which involve only the same 
lepton flavour £ for which we wish to compute the anomalous magnetic moment, and those which 
involve loops with leptons of different flavours, denoted collectively as £' [a = e 2 /4ir], 

a\ n „ ,„ „ /N (a 



«? ED = E^ (z) ! • (2.21) 



n>l n>2 



The second type of contribution, a^ ad , involves also quark loops. Their contribution is far from being 
limited to the short distance scales, and a^ ad is an intrinsically non perturbative quantity. From a 
theoretical point of view, this represents a serious difficulty. Finally, at some level of precision, the 
weak interactions can no longer be ignored, and contributions of virtual Higgs or massive gauge boson 
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degrees of freedom induce the third component aj . Of course, starting from the two loop level, a 
hadronic contribution to aj cak will also be present. The remainder of this presentation is devoted to 
a detailed discussion of these various contributions. 

Before starting this guided tour of the anomalous magnetic moments of the massive charged leptons 
of the standard model, it is useful to keep in mind a few simple considerations: 

• The anomalous magnetic moment is a dimensionless quantity. Therefore, the coefficients A n above 
are universal, i.e. they do not depend on the flavour of the lepton whose anomalous magnetic moment 
we wish to evaluate. 

• The contributions to ag of degrees of freedom corresponding to a typical scale M S> rri£ decouple 
[18], i.e. they are suppressed by powers of rrii/M. 7 

• The contributions to originating from light degrees of freedom, characterized by a typical scale 
m <C mi are enhanced by powers of \n{mi/m). At a given order, the logarithmic terms that do not 
vanish as m^/m — ► can often be computed from the knowledge of the lower order terms and of the 
(3 function through the renormalization group equations [24, 25, 26, 27]. 

These general properties already allow to draw several elementary conclusions. The electron being the 
lightest charged lepton, its anomalous magnetic moment is dominantly determined by the values of the 
coefficients A n . The first contribution of other degrees of freedom comes from graphs involving, say, at 
least one muon loop, which occurs first at the two-loop level, and is of the order of (m e /m ll ) 2 (a/Tr) 2 ~ 
10 -10 . The hadronic effects, i.e. "quark and gluon loops", characterized by a scale of ~ 1 GeV, or 
effects of degrees of freedom beyond the standard model, which may appear at some high scale M, 
will be felt more strongly, by a considerable factor (m M /m e ) 2 ~ 40 000, in a M than in a e . Thus, a e is 
well suited for testing the validity of QED at higher orders, whereas a M is more appropriate for testing 
the weak sector of the standard model, one of the main motivations for the BNL experiment, and 
possibly for detecting new physics. 



Exercises for section 2 

Exercise 2.1 

Show that the expression of the matrix element of the electromagnetic current given by Eqs. (2.4) and 
(2.5) indeed follows from the conditions stated. Show that the form factors Fi(k 2 ) in these equations 
are real if the current J p is hermitian. Work out the transformation properties of the different form 
factors under the operations of parity and time reversal. How many additional form factors are needed 
in order to describe the same matrix element if the assumption concerning the conservation of the 
current is dropped? 

Exercise 2.2 

Show that the current J p defined by Eq. (2.8) is conserved. 

7 In the presence of the weak interactions, this statement has to be reconsidered. Indeed, the necessity for the 
cancellation of the SU(2) x (7(1) gauge anomalies [19, 20, 21] transforms the decoupling of, say, a single heavy fermion 
in a given generation, into a somewhat subtle issue [22, 23], the resulting lagrangian being no longer renormalizable. 
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Exercise 2.3 

Find the term one needs to add to £j nt , and thus to J p , such as to generate a constant but nonzero 
form factor F4 at tree level. Show that in the non relativistic limit it induces an interaction term of 
the form F±(Qi)<t ■ (V A B) in a non uniform magnetic field. 

Exercise 2.4 

Show that the quantities A?(p',p) defined in Eq. (2.15) indeed project on the corresponding form 
factors Fi(k 2 ) through Eq. (2.14). Derive Eqs. (2.16) and (2.17). 

Exercise 2.5 

Work out the expression of the electromagnetic current in the standard model. 
Exercise 2.6 

Give the most general decomposition of the matrix element in Eq. (2.5) in the case of a massive 
Majorana neutrino [hint: see Refs. [28, 29] for a rather complete treatment]. 



3 Brief overview of the experimental situation 

3.1 Measurements of the magnetic moment of the electron 

The first indication that the gyromagnetic factor of the electron is different from the value g e = 2 
predicted by the Dirac theory came from the precision measurement of hyperfine splitting in hydrogen 
and deuterium [30] . The first measurement of the gyromagnetic factor of free electrons was performed 
in 1958 [31], with a precision of 3.6%. The situation began to improve with the introduction of 
experimental setups based on the Penning trap. Some of the successive values obtained over a period 
of forty years are shown in Table 1. Technical improvements, eventually allowing for the trapping of a 
single electron or positron, produced, in the course of time, an enormous increase in precision which, 
starting from a few percents, went through the ppm levels, before culminating at 4 ppb [parts per 
billion] in the last [32] of a series of experiments performed at the University of Washington in Seattle. 
The same experiment has also produced a measurement of the magnetic moment of the positron with 
the same accuracy, thus providing a test of CPT invariance at the level of 1CP 12 , 

9e-/9e+ = 1 + (0-5 ± 2.1) x 10- 12 . (3.1) 

Assuming invariance under CPT, the weighted average of the electron and positron anomalous mo- 
ments obtained in Ref. [32] gives [33, 34], 

a fP = 0.001159 652188 3(4 2). (3.2) 

An extensive survey of the literature and a detailed description of the various experimental aspects 
can be found in [35]. The earlier experiments are reviewed in [36]. 

3.2 Measurements of the magnetic moment of the muon 

The anomalous magnetic moment of the muon has also been the subject of quite a few experiments. 
The very short lifetime of the muon, r M = (2.19703 ± 0.00004) x 10~ 6 s, makes it necessary to proceed 
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Table 1: Some experimental determinations of the electron's anomalous magnetic moment a e with the 
corresponding relative precision. 



0.001 19(5) 


4.2% 


[37] 


0.001 165(11) 


1% 


[38] 


0.001 116(40) 


3.6% 


[31] 


0.001 160 9(2 4) 


2 100 ppm 


[39] 


0.001 159 622(27) 


23 ppm 


[40] 


0.001 159 660(300) 


258 ppm 


[41] 


0.001 159 657 7(3 5) 


3 ppm 


[42] 


0.001 159 652 41(20) 


172 ppb 


[43] 


0.001 159 652 188 4(4 3) 


4 ppb 


[32] 



in a completely different way in order to attain a high precision. The experiments conducted at CERN 
during the years 1968-1977 used a muon storage ring [for details, see [44] and references quoted therein]. 
The more recent experiments at the AGS in Brookhaven are based on the same concept. Pions are 
produced by sending a proton beam on a target. The pions subsequently decay into longitudinally 
polarized muons, which are captured inside a storage ring, where they follow a circular orbit in the 
presence of both a uniform magnetic field and a quadrupole electric field, the latter serving the purpose 
of stabilizing the orbits. The difference between the spin precession frequency and the orbit frequency 
is given by 

1 



u s - u c = <^ a„B 



a a + 



7 2 



/3AE , (3.3) 



where (3 is the velocity of the muons, and 7 is the corresponding Lorentz boost factor. Therefore, 
if 7 is tuned to its "magic" value 7 = ^ 1 + 1 ja^ = 29.3, the measurement of lo s — u) c and of the 
magnetic field B allows to determine a^. The spin direction of the muon is determined by detecting 
the electrons or positrons produced in the decay of the muons with an energy greater than some 
threshold energy E t . The number of detected electrons N e (t) decreases exponentially with time, the 
time constant being set by the muon's lifetime 7T M c in the laboratory frame, and is modulated by the 
frequency uj s — uj c , 

N e (t) = iVo(£ t )e-*/ 77 >{l + A(E t ) cos[(u; s - to c )t + (j>(E t )}} . (3.4) 
The observation of this time dependence thus provides the required measurement of uj s — uj c - 



Several experimental results for the anomalous magnetic moment of the positively charged muon, 
obtained at the CERN PS or, more recently, at the BNL AGS, are recorded in Table 2. Notice that 
the relative errors are measured in ppm units, to be contrasted with the ppb level of accuracy achieved 
in the electron case. The four last values in Table 2 were obtained by the E821 experiment at BNL. 
They show a remarkable stability and a steady increase in precision, and now completely dominate 
the world average value. Further data, for negatively charged muons 8 are presently being analysed. 

8 The CERN experiment had also measured a M - = 0.001 165 937(12) with a 10 ppm accuracy, giving the average value 
a M = 0.001 165 924(8.5), with an accuracy of 7 ppm. 
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Table 2: Determinations of the anomalous magnetic moment of the positively charged muon from the 
storage ring experiments conducted at the CERN PS and at the BNL AGS. 



0.001 16616(31) 


265 ppm 


[45] 


0.001 165 895(27) 


23 ppm 


[46] 


0.001 165 911(11) 


10 ppm 


[47] 


0.001 165 925(15) 


13 ppm 


[48] 


0.001 165 9191(5 9) 


5 ppm 


[49] 


0.001 165 920 2(16) 


1.3 ppm 


[50] 


0.001 165 920 3(8) 


0.7 ppm 


[51] 



The aim of the Brookhaven Muon (g - 2) Collaboration is to reach a precision of 0.35 ppm, but this 
will depend on whether the experiment will receive financial support to collect more data or not 9 . 

For completeness, one should mention that Eq. (3.3) is only correct as long as the muon has no electric 
dipole moment. If this is not the case, the more general relation, 

u s - u c = — \ 

which holds for (3 ■ E = (3 ■ B = 0, has to be used. The additional term proportional to induces an 
oscillation of the muon spin with respect to the plane of motion 10 . In the standard model, given the 
experimental precision and the intensities of the fields used in the experiment, it is quite legitimate to 
use the formula (3.3) instead of (3.5). However, in case a discrepancy arises between the experimental 
value of a M and the standard model prediction, the difference could be induced by non standard 
contributions to either the anomalous magnetic moment or the electric dipole moment d^, see the 
discussion in [54]. 



a u + 



1-7 2 



2d 

/3AE[ - -^{/3AB + E} 



(3.5) 



3.3 Experimental bounds on the anomalous magnetic moment of the r lepton 

As already mentioned, the very short lifetime of the r precludes a measurement of its anomalous 
magnetic moment following any of the techniques described above. Indirect access to a T is provided 
by the reaction e + e~ — > t + t~j. The results obtained by OPAL [56] and L3 [57] at LEP only lead to 
very loose bounds, 

-0.052 < a T < 0.058 (95%C.L.) 

-0.068 <a T < 0.065 (95%C.L.) , (3.6) 
9 At the time of writing, the prospects in this respect are unfortunately rather dim. 

10 A measurement of the electric dipole moment of the muon was actually performed by the CERN experiment [47], 
with the result d M = (3.7 ± 3.4) x 1CP 19 ecm. An even smaller value, d M < 9.1 x 1CP 25 ecm, can be infered from the 
experimental value of the electric dipole moment of the electron [52, 53], d e = 1.8(1.2)(1.0) x lCF 27 ,ecm, and assuming 
a scaling law d M ~ 7^~d e . Such a scaling law holds within the standard model, but not in models with flavour violating 
interactions, see for instance [54]. For a proposal to measure d M at the level of ~ lCT 24 ,ecm, see [55]. 
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respectively. 



We shall now turn towards theory, in order to see how the standard model predictions compare with 
these experimental values. Only the cases of the electron and of the muon will be treated in some 
detail. The theoretical aspects as far as the anomalous magnetic moment of the r are concerned are 
discussed in [58] and in the references quoted therein. 



4 The anomalous magnetic moment of the electron 

We start with the anomalous magnetic moment of the lightest charged lepton, the electron. Since 
the electron mass m e is much smaller than any other mass scale present in the standard model, the 
mass independent part of o^ ED dominates its value. As mentioned before, non vanishing contributions 
appear at the level of the loop diagrams shown in Fig. 1 . 

H 

QED + y<^^\^ + y(^J^ + 

Figure 1: The perturbative expansion of T p (p' ,p) in single flavour QED. The tree graph gives F\ = 1, 
F2 = 0, whereas F^{k 2 ) and F^{k 2 ) vanish identically to all orders in pure QED. The one loop vertex 
correction graph gives the coefficient A\ in Eq. (2.21). The cross denotes the insertion of the external 
field. 



4.1 The lowest order contribution 

The one loop diagram gives 



rvalue* 



He) 5 



A 

(27T)*' 
i) 



f'+fi- 



— r 



(4.1) 



The photon propagator has been written in a Lorentz type gauge, corresponding to a covariant gauge 
fixing term —{d ■ A) 2 /2£. Let us for a moment concentate on the £ - dependence of this one loop 
expression. Recall that T p (p',p) has eventually to be inserted between the spinors u(p') and u(p). 
Then, the gauge dependence of the integrand is given by 



(i-Ou(pO^-tt 



m,. 



Mp) = (1 - £)u(p')«7/>Hi(p) 



'+ $ — m e ' p + 

and thus affects F\{k 2 ), but not F2(k 2 ). For evaluating the latter, one may therefore take, say, £ = 1 
for convenience. The form factor F2(k 2 ) is then obtained by using Eqs. (2.14) and (2.15) and, upon 
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evaluating the corresponding trace of Dirac matrices, one finds 



F 2 (k 2 



32m{ 



1 1 loop 



le 



d A q 



1 



k 2 (k 2 - Am 2 e ) 2 J ( 27r ) 4 {p' + q) 2 -m 2 (p + q) 2 



m 



1 

2 ^2 



-3k 2 (p ■ q) 2 + 2k 2 m 2 (p ■ q) + k 2 m 2 e q 2 - m 2 e (k ■ q) 



(4.2) 



Then follow the usual steps of introducing two Feynman parameters, of performing a trivial change 
of variables and a symmetric integration over the loop momentum q, so that one arrives at 



64m 2 



(k 2 - Am 2 ) 2 



(2tt) 4 (q 2 - Tl 2 f 



2x(l — x)m 



x y (k ) + m e k x 3xy — y + —x 



2m 2 



ir 2 {k 2 - Am 2 ) 2 

x 



L dxx L dv h 



2x{\ — x)m\ — -x 2 y 2 (k 2 ) 2 + m 2 e k 2 x ( 3xy 



with 



n 2 = x 2 y{\ - y){2m 2 e - k 2 ) + x 2 y 2 m 2 e + x 2 (l - yfm\ 
As expected, the limit k 2 —>■ can be taken without problem, and gives 

1 a 



a e \i loop = F 2 (0)\ 1 loop 



2 7T 



(4.3) 
(4.4) 

(4.5) 



Let us stress that although the integral (4.1) diverges, we have obtained a finite result for ^(A; 2 ), and 
hence for a e , without introducing any regularization. This is of course expected, since a divergence in 
^2(0) would require that a counterterm of the form given by the second term in A n t, see Eq. (2.8), 
be introduced. This would in turn spoil the renormalizability of the theory. In fact, as is well known, 
the divergence lies in i'i(O), and goes into the renormalization of the charge of the electron. 




+ sym 

Figure 2: The Feynman diagrams which contribute to the coefficient A 2 in Eq. (2.21). 
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4.2 Higher order mass independent corrections 

The previous calculation is rather straightforward and amounts to the result 

Ai = \ (4.6) 

first obtained by Schwinger [59]. Schwinger's calculation was soon followed by a computation of 
A 2 [60], which requires the evaluation of 7 graphs, representing five distinct topologies, and shown in 
Fig. 2. Historically, the result of Ref. [60] was important, because it provided the first explicit example 
of the realization of the renormalization program of QED at two loops. However, the value for A 2 was 
not given correctly. The correct expression of the second order mass independent contribution was 
derived in [61, 62, 63] (see also [64, 65]) and reads 11 

197/1 \ 3 

A * = l44 + (2- 31n2 ) C(2) + 4 C(3) 

= -0.328478 965... (4.7) 

oo 

with C{p) = V nP ' C(2) = ^ 2 /6- The occurence of transcendental numbers like ((p) is a general 

n=l 

feature of higher order calculations in perturbative quantum field theory. The pattern of these tran- 
scendental in perturbation theory and the structure of the renormalization algorithm have also been 
put in relationship with other mathematical structures, like knot theory and braids [66], Hopf algebras 
[67] and non commutative geometry [68]. 

The analytic evaluation of the three loop mass independent contribution to the anomalous magnetic 
moment required quite some time, and is mainly due to the dedication of E. Remiddi and his cowork- 
ers during the period 1969-1996. There are now 72 diagrams to consider, involving many different 
topologies, see Fig. 3. 



The calculation was completed [69] in 1996, with the analytical evaluation of a last class of diagrams, 
the non planar "triple cross" topologies. The result reads 12 



83 



215 



100 



A 3 = _7r'C(3) - — C(5) + — 



72 



239 



:7T* + 



24 
139 



«4 + ^ln 4 2 



298 



2160 
1.181241456... 



C(3) - — vr 2 In 2 + 
18 v ' 9 



17101 
810 



-7r 2 ln 2 2 
24 

28259 



5184 



(4.8) 



where 13 a p = ^ ^ n p - The numerical value extracted from the exact analytical expression given 



n=l 



above can be improved to any desired order of precision. 

11 Actually, the experimental result of Ref. [38] disagreed with the value A2 = —2.973 obtained in [60], and prompted 
theoreticians to reconsider the calculation. The result obtained by the authors of Refs. [61, 62, 63] reconciled theory 
with experiment. 

12 The completion of this three-loop program can be followed through Refs. [70]- [75] and [69]. A description of the 
technical aspects related to this work and an account of its status up to 1990, with references to the corresponding 
literature, are given in Ref. [76]. 

13 The first three values are known to be 01 = In 2, a 2 = Li 2 (l/2) = (£(2) - In 2 2)/2, a 3 = |C(3) - s<(2) In 2 + § In 3 2 
[76]. 
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12 4 6 

Figure 3: The 72 Feynman diagrams which make up the coefficient ^3 in Eq. (2.21). 



In parallel to these analytical calculations, numerical methods for the evaluation of the higher order 
contributions were also developed, in particular by Kinoshita and his collaborators (for details, see 
[77]). The numerical evaluation of the full set of three loop diagrams was achieved in several steps [78]- 
[84]. The value quoted in [84] is A3 = 1.195(26), where the error comes from the numerical procedure. 
In comparison, let us quote the value [85, 77] A3 = 1.176 11 (42) obtained if only a subset of 21 three 
loop diagrams out of the original set of 72 is evaluated numerically, relying on the analytical results 
for the remaining 51 ones, and recall the value ^3 = 1.181241456... obtained from the full analytical 
evaluation. The error induced on a e due to the numerical uncertainty in the second, more accurate, 
value is still A(a e ) = 5.3 x 10 -12 , whereas the experimental error is only A(a e )| 0X p = 4.3 x 10 -12 . 
This discussion shows that the analytical evaluations of higher loop contributions to the anomalous 
magnetic moment of the electron have a strong practical interest as far as the precision of the theoretical 
prediction is concerned, and which goes well beyond the mere intellectual satisfaction and technical 
skills involved in these calculations. 14 

At the four loop level, there are 891 diagrams to consider. Clearly, only a few of them have been 
evaluated analytically [86, 87]. The complete numerical evaluation of the whole set gave [85] A4 = 
— 1.434(138). The developement of computers allowed subsequent reanalyses to be more accurate, 
i.e. A4 = —1.557(70) [88]. Until recently, the "latest of [these] constantly improving values" was 
[5] A/± = —1.509 8(38 4). This calculation certainly represents a formidable task, and requires many 
elaborate technical tools. A descriptive account can be found in [77]. Let us mention, for completeness, 
that efforts to improve upon the evaluation of A4 are presently being pursued. Thus, a mistake has 
recently been found in an earlier computer code used for the evaluation of a subclass of four loop 
diagrams [89], whose contribution to A4 was A^-jvCd) = —0.7503(60) [for the precise meaning of 
the notation, we refer the reader to Refs. [77] and [85]]. The corrected value reads instead [89] 

14 It is only fair to point out that the numerical values that are quoted here correspond to those given in the original 
references. It is to be expected that they would certainly improve if today's numerical possibilities were used. 



14 



^4;iv(d) = —0.99072(10) [note also the impressive improvement in the precision]. By itself, this 
correction induces a 16% downward shift of the value of A4, a far from trivial modification [see below]. 
At the present stage, the value of A4 reads 

A A = -1.750 2(384). (4.9) 



Needless to say, so far the five loop contribution A5 is unknown territory. On the other hand, (a/V) 5 ~ 
7 x 10~ 14 , so that one may reasonably expect, in view of the present experimental situation, that its 
knowledge is not yet required. 

4.3 Mass dependent QED corrections 

We now turn to the QED loop contributions to the electron's anomalous magnetic moment involving 
the heavier leptons, /j, and r. The lowest order contribution of this type occurs at the two loop level, 
0(a 2 ), and corresponds to a heavy lepton vacuum polarization insertion in the one loop vertex graph, 
cf. Fig. 4. Quite generally, the contribution to ag arising from the insertion, into the one loop vertex 
correction, of a vacuum polarization graph due to a loop of lepton £', reads [90, 91] 15 



B 2 



re oi\ 1 f j /1 4mf,s + 2m2, 1 x 2 (l - x) 

(1,1') = - / ds\ 1 * =-£- / dx ; f-j-. (4.10 

3 Jim? V 8 S 2 Jo X 2 + (1 - x) ^ 



If miii 3> mi, the second integrand can be approximated by x m 2 /s, and one obtains [93] 

B 2 (£,£') = f— V + •■• , m t >^>m e . (4.11) 



45 \m£ 

The complete expansion of B 2 (£,£') for m^i S> ni£ can be found in [71], from which we quote 
B 2 (£,£>) = If^V + lf^Vlnf^U 9 (m ^' 



45 \me> J 70 \me> J \mg> J 19600 
+ f^\ 6 ln (™>t \ 131 fm^" 



315 \m(i ) \mf i J 99225 \m£i 
+ l^Yln fme 



mil J \ mi 



(4.12) 



Numerically, this translates into [the values for the masses that were used read m e = 0.510 998 902(21) 
MeV, m T = 1776.99^ [34], and m^/m e = 206.768 277(24) [94]] 

B 2 (e,n) = 5.197 x 10~ 7 

S 2 (e,r) = 1.838 xl0~ 9 . (4.13) 



15 A trivial change of variable on s brings the expression (4.10) into the form given in [90, 91]. Furthermore, the 
analytical result obtained upon performing the double integration is available in [92]. 
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Figure 4: The insertion of a muon vacuum polarization loop into the electron vertex correction (left) 
or of an electron vacuum polarization loop into the muon vertex correction (right). 



For later use, it is interesting to briefly discuss the structure of Eq. (4.10). The quantity which appears 
under the integral is related to the cross section for the scattering of a pair into a pair 
at lowest order in QED, 



QED 



so that 



where 16 



Ana 2 , 

l r°° ds 

3 JAm 2 tl S 



1- ^(s + 2m|,). 



K{s)R^'\s) 



K{s) 



1 x 2 (l-x) 

o x 2 + (1 - x) A 

v ' mi 

£+ ■£--+(£• ')+(£')- 



and R^'\s) is the lowest order QED cross section o-q ED 
of the cross section of the reaction e + e~ — > [i + for s 3> m 2 cr, 



(4.14) 
(4.15) 

(4.16) 

(s) divided by the asymptotic form 



' oo 



(s) 



Ana'' 
3s 



The three loop contributions with different lepton flavours in the loops are also known analytically 
[95, 96]. It is convenient to distinguish three classes of diagrams. The first group contains all the 
diagrams with one or two vacuum polarization insertion involving the same lepton, [i ot r, of the type 
shown in Fig. 5. The second group consists of the leptonic light-by-light scattering insertion diagrams, 
Fig. 6. Finally, since there are three flavours of massive leptons in the standard model, one has also 
the possibility of having graphs with two heavy lepton vacuum polarization insertions, one made of a 
muon loop, the other of a r loop. This gives 



B 3 (e, 



Bt P -\e;») + Br>(e;r) + B^^A*) + B^"^(e;r) + ^(ej/i.r). (4.17) 



>( v -p-). 



? (LxL) 



(LxL), 



?( V -P-)/ 



The analytical expression for B^' p '^ (e; fj) can be found in Ref. [95] , whereas [96] gives the correspond- 
ing result for B^ xL \e;fi). For practical purposes, it is both sufficient and more convenient to use 
their expansions in powers of m e /mn, 



B 



(v.p.), 



e\fi) 



m 



135 V Tn e 



10117 



45 77 24300 



16 Explicit expressions for K(s) are also available, but for many purposes, the integral representation given here turns 
out to be more convenient. 
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+ 



J9_ 2 K 
2520 V m e 



14233 /m„ 
132300 V m P 



+ 



49 
768 



C(3)- 



11 

945 



7T^ + 



2976691 



296352000 



+ O 



-0.000 021768... 



(4.18) 




Figure 5: Three loop QED corrections with insertion of a heavy lepton vacuum polarization which 



make up the coefficient B% p '\e; fi). 



and 



B 



(LxL) 



(e;/i) 



2 r 



771. , 



5«s) - ^ 

2 SV ; 16 



+ 



nif. 



777, 



161 2 / 77^ 
810 Ue 



+ 



6' 



0.0000143945... 



16189 /77V 
48600 V 77i P 



161 
9720' 



-7T 



831931 
972000 



(4.19) 



The expressions for ZJg (e; r) and B^ x (e;r) follow upon replacing the muon mass by m T . 
This again gives a suppression factor (m^/mr) 2 , which makes these contributions negligible at the 
present level of precision. For the same reason, B^' p '\e; (j,, r) can also be discarded. 

4.4 Other contributions to a e 

In order to make the discussion of the standard model contributions to a e complete, there remains 
to mention the hadronic and weak components, cig ad and a™ eak , respectively. Their features will be 
discussed in detail below, in the context of the anomalous magnetic moment of the muon. I therefore 
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/ 

6 

Figure 6: The three loop QED correction with the insertion of a heavy lepton light-by-light scattering 
subgraph, corresponding to the coefficient B^ xL \e; fi). 

only quote the numerical values 1 

a^ d = 1.67(3) x 1(T 12 , (4.20) 

and [97] 

a » eak = 0.030 x 10~ 12 . (4.21) 



4.5 Comparison with experiment and determination of a 

Summing up the various contributions discussed so far gives the standard model prediction [4, 5, 8, 

2 /^,\3 / „, \ 4 



a e SM 



0.5- - 0.32847844400 ( -\ + 1.181 234 017 ( -\ - 1.750 2(384) ( -\ + 1.70(3) x 10~ 12 . 

(4.22) 

In order to obtain a number that can be compared to the experimental result, a sufficiently accurate 
determination of the fine structure constant a is required. The best available measurement of the 
latter comes from the quantum Hall effect [33], 

a - l (qH) = 137.036 003 00(2 70) (4.23) 

and leads to 

af M (qH) = 0.001159 652146 5(240), (4.24) 

about six times less accurate than the experimental value of Eq. (3.2), a| xp = 0.001 159 652 188 3(4 2). 
On the other hand, if one excludes other contributions to a e than those from the standard model 
considered so far, i.e. if one identifies a^ M with a| xp , then the value of a| xp as given in Eq. (3.2) 
provides the best determination of a to date [89, 8], 

a _1 (a e ) = 137.035 998 75(52), (4.25) 

at least to the extend that one may reasonably believe that all theoretical errors are under control. 
Now, as we have seen earlier, the value of has recently been corrected [89]. The analysis presented 

17 1 reproduce here the values given in [4, 5], except for the fact that I have taken into account the changes in the 
value of the hadronic light- by- light contribution to a M , see below, for which I take a^ LxL ' = +8(4) x 10~ 10 , and which 
translates into ai LxL ' ~ a{ t LxL ' (m e /m M ) 2 = 0.02 x 10~ 12 . 
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here incorporates the changes brought forward by the analysis of Ref. [89]. It lowers the prediction 
for a e by ~ 7.0 x 1CT 12 if one uses the value (4.23) of a, or equivalently reduces the value of a~ l (a e ) 
by one and a half standard deviation. It is very likely that the completion of the analysis begun in 
Ref. [89] will lead to a more accurate determination of A4, and further changes in the numbers quoted 
here are to be expected. 



Exercises for section 4 

Exercise 4-1 

Reproduce the steps that lead from Eq. (4.1) to Eq. (4.5). 



5 The anomalous magnetic moment of the muon 

In this section, the theoretical aspects concerning the anomalous magnetic moment of the muon are 
discussed. Since the muon is much heavier than the electron, a„ will be more sensitive to higher 
mass scales. In particular, it is a better probe for possible degrees of freedom beyond the standard 
model, like supersymmetry. The drawback, however, is that will also be more sensitive to the non 
perturbative strong interaction dynamics at the ~ 1 GeV scale. 

5.1 QED contributions to a M 

As already mentioned before, the mass independent QED contributions to are described by the 
same coefficients A n as in the case of the electron. We therefore need only to discuss the coefficients 
B n fa£') associated with the mass dependent corrections. 



For m e > <C m e , Eq. (4.10) gives [90, 91, 92, 71] 

• (5-1) 



w,n = ^ m ( - 1 + 4 — - 4 f — ) 2 ln f — ) + 3 f — ) 2 + ° 

3 \m v ) 36 4 m e \m e J \m v ) \m e J 



The complete expansion in powers of m^jm^ can again be found in [71]. In the case of B2(/j,;t), 
one may use the expression given in Eq. (4.12). Upon using the values m e = 0.510 998 902(21) MeV, 
rrif, = 105.658 357(5) MeV, m T = 1776.99±g;|§ MeV [34], and m M /m e = 206.768 277(24) [94] the 
corresponding numbers read 

B 2 (fi;e) = 1.094 258 300(38) (5.2) 

B 2 (h;t) = 0.000 078 064(25). (5.3) 

Although these numbers follow from an analytical expression, there are uncertainties attached to them, 
induced by those on the corresponding values of the ratios of the lepton masses. 

The three loop QED corrections decompose as 

B 3 (M) = i#" P V;e) + B$ v - p V;r) + B^ xL) fa e) + B^far) + B^ fa e,r) , (5.4) 
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with [95, 96] 
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(5.6) 



while B^' p '\fi; r) and B^ xL \fi;r) are derived from £?3 V ' p '^(e; /x) and from B^ xL \e; fi) , respectively, 
by trivial substitutions of the masses. Furthermore, the graphs with mixed vacuum polarization 
insertions, one electron loop, and one r loop, are evaluated numerically using a dispersive integral 
[71, 95, 98]. 
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Numerically, one obtains 

B^ p -\n;e) = 1.9204551(2) 

B^ xL \fi;e) = 20.947924 7(7) 

^3 V ' P ' ) (m;t) = -0.001782 3(5) 

B^ xL \v,t) = 0.002142 9(7) 

B£*-\fi;e,T) = 0.000 527 7(2). (5.7) 

Notice the large value of B^ xIj \fi; e), due to the occurence of terms involving factors like ln(m^/m e ) ~ 
5 and powers of tt. Such a large contribution, first obtained numerically in Ref. [78], allowed to explain 
a discrepancy of 1.7 a between the theoretical value and the experimental measurement of Ref. [45]. 
Finally, several pieces of B^' p '\fi; e) had already been worked out earlier, in Refs. [24, 93, 99, 100, 101] 

The contributions at fourth order in a have been obtained numerically in Ref. [102]. They must be 
corrected for the change in A4 obtained in [89]. At the next order, no full calculation, even through 
numerical techniques, is available. Specific contributions, for instance those enhanced by powers of 
ln(m e /m M ) times powers of n, have been evaluated [103, 104, 105, 106]. 

Putting all these contributions together leads to the expression 

a Q ED = 0.5^ + 0.765 857399(45) Q + 24.050 509 5(2 3) + 125.08(41) Q +930(170) 

(5*8) 

Upon inserting the value of a obtained from the anomalous magnetic moment of the electron in Eq. 
(4.25), one finds 

a QED = n 6 5 8 470.35(28) x 10~ 10 . (5.9) 
5.2 Hadronic contributions to a M 

On the level of Feynman diagrams, hadronic contributions arise through loops of virtual quarks and 
gluons. These loops also involve the soft scales, and therefore cannot be computed reliably in per- 
turbative QCD. We shall decompose the hadronic contributions into three subsets: hadronic vacuum 
polarization insertions at order a 2 , at order a 3 , and hadronic light-by-light scattering, 

a£ ad = of - v - p - 1} + of - v - p - 2) + of" LxL) (5.10) 
5.2.1 Hadronic vacuum polarization 

We first discuss af ' v ' p ' which arises at order 0(a 2 ) from the insertion of a single hadronic vac- 
uum polarization into the lowest order vertex correction graph, see Fig. 7. The importance of this 
contribution to a M is known since long time [107, 108]. 



There is a very convenient dispersive representation of this diagram, similar to Eq. (4.10) 



a (h.v.p. 1) = 4a 2 / ^£ K ,s £ Imn(s) 
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Figure 7: The insertion of the hadronic vacuum polarization into the one loop vertex correction, 
corresponding to a^ h ' v ' p ' . 

= \(-) 2 r 2 -K(s)R^(s), (5.11) 

Here, n(s) denotes the hadronic vacuum polarization function, defined as 18 

(q»qu ~ q\uMQ 2 ) = i J dtxe*** {n\T{j^x)j u (0)}\Q) , (5.12) 

with j p the hadronic component of the electromagnetic current, Q 2 = —q 2 > for q^ spacelike, and 
the QCD vacuum. The function K(s) was defined in Eq. (4.16), and R had (s) stands now [see 
however below] for the cross section o~q +£ ^ had (s) of e + e~ — ► hadrons, at lowest order in a, divided 
by °1» e ~^ +At (s) = 4 g° . A first principle computation of this strong interaction contribution is far 
beyond our present abilities to deal with the non perturbative aspects of confining gauge theories. This 
last relation is however very interesting because it expresses a^ h ' v ' p ' ^ through a quantity that can be 
measured experimentally. In this respect, two important properties of the function K(s) deserve to 
be mentioned. First, it appears from the integral representation (4.16) that K{s) is positive definite. 
Since R had is also positive, one deduces that a^ h v ' p ' itself is positive. Second, the function K(s) 
ill 



decreases as m 2 /3s as s grows, so that it is indeed the low energy region which dominates the integral. 



Explicit evaluation of a^ h ' v p ' using available data actually reveals that more than 80% of its value 
comes from energies below 1.4 GeV. This observation is rather welcome, since the spectral density 
ImII(s) can also be extracted from data on the hadronic decays of the r lepton in this energy region. 

Finally, the values obtained in this way for a^ h v ' p ' 1 ' have evolved in time, as shown in Table 3. This 
evolution is mainly driven by the availability of more data, and is still going on, as the last entries of 
Table 3 show. In order to match the precision reached by the latest experimental measurement of a^, 
^(h.v.p. l) neec j s £ k e kno^n a i \% Besides the very recent high quality e + e~ data obtained by 
the BES Collaboration [109] in the region between 2 to 5 GeV, and by the CMD-2 collaboration [110] 
in the region dominated by the p resonance, the latest analyses sometimes also include or use, in the 
low-energy region, data obtained from hadronic decays of the r by ALEPH [111], and, more recently, 
by OPAL [112] and CLEO [113, 114]. We may notice from Table 3 that the precision obtained by using 

18 Actually, II(s) defined this way has an ultraviolet divergence, produced by the QCD short distance singularity of the 
chronological product of the two currents. However, it only affects the real part of II(s). A renormalized, finite quantity 
is obtained by a single subtraction, Il(s) — 11(0). 
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e+e" data alone has become comparable to the one achieved upon including the r data. However, 
one of the latest analyses reveals a troubling discrepancy between the purely e+e" and the r based 
evaluations. Additional work is certainly needed in order to resolve these problems 19 . Further data 
are also expected in the future, from the KLOE experiment at the DAPHNE e+e~ machine [124], or 
from the B factories BaBar [125] and Belle. For additional comparative discussions and details of the 
various analyses, we refer the reader to the literature quoted in Table 3. Averaging the results from 
the three most recent analysis with e+e" data only, from Refs. [121], [122], and [123], gives [8] 

a (h.v. P . i)( e + e -) = 68 3 8 (75) x iQ- 11 . (5.13) 



Table 3: Some of the recent evaluations of ' v p ' ' x 10 11 from e+e~ and/or r-decay data. Note that 
the authors of Ref. [120] also use space like data for the pion form factor F n (t); furthermore, they use 
a preliminary version of the CMD-2 data, that is now superseded by [110]. 



7024(153) 


[115] 


e+e - 


7026(160) 


[116] 


e+e - 


6950(150) 


[117] 


e+e - 


7011(94) 


[117] 


r, e + e", 


6951(75) 


[118] 


r, e+e", QCD 


6924(62) 


[119] 


r, e+e", QCD 


7016(119) 


[58] 


e+e", QCD 


7036(76) 


[58] 


r, e+e - , QCD 


7002(73) 


[120] 


e+e - , incl. BES-II data, F n (t) 


6836(86) 


[121] 


e+e", incl. BES-II and CMD-2 data 


6847(70) 


[122] 


e+e", incl. BES-II and CMD-2 data 


7090(59) 


[122] 


r, e+e", incl. BES-II data 


6831(62) 


[123] 


e+e", incl. BES-II and CMD-2 data 



Let us briefly mention here that it is quite easy to estimate the order of magnitude of a^ h ' v ' p ' ^. For 
this purpose, it is convenient to introduce still another representation [2, 126], which relates a^ h ' v p ' ^ 
to the hadronic Adler function A(Q 2 ), defined as 20 



, 2 dH(Q 2 ) f°° Q 2 1 



by 



Am - -Q - j o it - I mI I( t ) , (5.14) 

= 2 , 2 Q^'| (1 _ x)(2 _ lM ^ m ^. (5 , 5) 



19 For a possible explanation, see 



20 Unlike 11(f) itself, A(Q 2 ) if free from ultraviolet divergences. 
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A simple representation of the hadronic Adler function can be obtained if one assumes that Imil(i) 
is given by a single, zero width, vector meson pole, and, above a certain threshold so, by the QCD 
perturbative continuum contribution, 

i ImII(t) = | flMlbit - M 2 V ) + | ^ [1 + 0{a s )\ 6(t - s ) (5.16) 

The justification [127] for this type of minimal hadronic ansatz can be found within the framework 
of the large-iVc limit [128, 129] of QCD, see Ref. [127] for a general discussion and a detailed study 
of this representation of the Adler function. The threshold sq for the onset of the continuum can be 
fixed from the property that in QCD there is no contribution in 1/Q 2 in the short distance expansion 
of A(Q 2 ), which requires [127] 

2tiM 2 =^s (l + l^ + O(a 2 )). (5.17) 

This then gives a^' v ' p ' ^ ~ (570 ± 170) x 10~ 10 , which compares well with the more elaborate data 
based evaluations in Table 3, even though this simple estimate cannot claim to provide the required 
accuracy of about 1%. 




Figure 8: Higher order corrections containing the hadronic vacuum polarization contribution, corre- 
sponding to a^ h ' v p ' 2 ^ . 

We now come to the C(a 3 ) corrections involving hadronic vacuum polarization subgraphs. Besides 
the contributions shown in Fig. 8, another one is obtained upon inserting a lepton loop in one of the 
two photon lines of the graph shown in Fig. 7. Taken all together, these can again be expressed in 
terms of i? had [130, 3, 98] 

a (h.v. P . 2) = 1 ("A 3 f°° ^(2) (s)jR had (s) _ (518) 

3 \TT J J4M% S 

The value obtained for this quantity is [98] 

a (h.v.p. 2) x 10 11 = _ 1Q1 ± g _ ( 5 19 ) 

The expression of KW(s) is given in the references quoted above. 
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Figure 9: A higher order correction containing the hadronic vacuum polarization contribution, and 



which is included in a^ h ' v p ' 1 \ 



There is actually another C(a 3 ) correction, namely the one obtained upon attaching a virtual photon 
line with both ends to the hadronic blob in Fig. 7, see Fig. 9. On the other hand, a^ h v ' p ' involves in 
principle data corrected for all electromagnetic effects. Whereas radiative corrections in the leptonic 
initial state and vacuum polarization effects in the photon propagator certainly can be accounted for, 
there is at present no way to handle in a model independent way electromagnetic corrections in the 
hadronic final state. These, on the other hand, contribute, together with final states containing an 
additional photon, to the O(o 3 ) contribution we have just been mentioning. It has become customary 
to include all these effects into a|j h v ' p ' 1 \ where it is then to be understood that R had (s) in Eq. (5.11) 
actually stands for 

^had (s) = ^0 + ^2 W+V W (5 20) 

where 0Q +e ^ had (s) + cr| +e ^ had (s) denotes the cross section for e + e~ — ► had beyond leading order 
in the expansion in powers of the fine structure constant a. The values given in Table 3 correspond 
to the definition (5.20). It also should be stressed that the next-to-leading cross section af e ^ had (s), 
as well as the radiative cross section o"q e ~* had+7 (s), are infrared divergent quantities, and only their 
sum is actually well defined. 

5.2.2 Hadronic light-by-light scattering 

We now discuss the so called hadronic light-by-light scattering graphs of Fig. 10. 

The contribution to T p (p',p) of relevance here is the matrix element, at lowest nonvanishing order in 
the fine structure constant a, of the light quark electromagnetic current 

2 1 - 1 

j P ( x ) = ^(u-f p u)(x) - -{d-y p d){x) - -(s-/ p s)(x) (5.21) 

between fj,~ states, 

(-ie)u(p')T p h - L * L Xp',p)u(p) ee </rV)l(ie)j,(o)|/r(p)) 

Ai r d% h) 3 



(27r)4y (2tt)4 qfqUqi+12-ky 
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permutations 



Figure 10: The hadronic light-by-light scattering graphs contributing to aj^' LxL \ 



(p' — Qi) 2 — m2 (p 1 — Qi — Q2) 2 — fn 2 
x(-iefu(p')^{p"- 4x + m)Y(p"- ix- d 2 + m)j x u(p) 
x (ie) 4 n AtJyAp (Q , i, <?2, k- qi-q 2 ) , (5.22) 



with ky, = {p' — p) M and 

IW P (<7i,<72,g3) 



x { n 1 T{j^ Xl )j u (x 2 )j x (x 3 )j p (o)} 1 n ) 



(5.23) 



the fourth-rank light quark hadronic vacuum-polarization tensor, | £1 ) denoting the QCD vacuum. 
Since the flavour diagonal current j^ix) is conserved, the tensor Tl piU \ p {qi, q 2 , q 3 ) satisfies the Ward 
identities 

{qi;q2;qt,(Qi + Q2 + q3) p }n pi/Xp ( qi ,q 2 ,q 3 ) = {0;0;0;0}. 

This entails that 21 



Q(p / )r (h. LxL )(p / jp)u(p) = fi(p ') 



7^ 



(h. LxL) (A;2) + J_ afrrk rp^ LxL) (fc2; 



u(p) 



(5.24) 
(5.25) 



as well as F p h ' LxL \p' ,p) = k T T p h T LxL \p' ,p) with 



u(p')r£ LxL V,p)u(p) = -ie e 



(f 4 gi /" (i 4 g 2 



1 



(2tt)4 7 (2vr)4 ^Ife+^-fc) 2 
1 1 



(p' — Qi) 2 ~ m2 ip' ~ Qi ~ Q2) 2 — m2 
xu(p')^(p"- 4i + m)Y{p"- i\- d 2 + m)j x vi{p) 
d 



(5.26) 



Following Ref. [78] and using the property k p k a \x{p')T p h a LxL \p' ,p)u(p) = 0, one deduces that 

F (h. LXL) 

(0) = and that the hadronic light-by-light contribution to the muon anomalous magnetic 

21 We use the following conventions for Dirac's 7-matrices: {7^,7,,} = 2r/ M „, with rj^ the flat Minkowski space metric 

of signature (H ), a^ v = (i/2)[y^, 7,,], 75 = i7°7 1 7 2 7 3 , whereas the totally antisymmetric tensor e^ vp a is chosen 

such that £oi23 = +1. 
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moment is equal to 

a (h. LXL) _ ^(h. LXL) (0) = _l_ tr | ( ^ + m)[7 p )r]( ^ + m)r (h. LxL )(pjp) j (537) 

This is about all we can say about the QCD four-point function H^XpiQi, Q2, $3)- Unlike the hadronic 
vacuum polarization function, there is no experimental data which would allow for an evaluation of 
Q (h. LxL)^ rp^ e ex j s ^j n g estimates regarding this quantity therefore rely on specific models in order to 
account for the non perturbative QCD aspects. A few particular contributions can be identified, see 
Fig. 11. For instance, there is a contribution where the four photon lines are attached to a closed 
loop of charged mesons. The case of the charged pion loop with pointlike couplings is actually finite 
and contributes ~ 4 x 10 -10 to [131]. If the coupling of charged pions to photons is modified 
by taking into account the effects of resonances like the p, this contribution is reduced by a factor 
varying between 3 [131, 132] and 10 [133], depending on the resonance model used. Another class of 
contributions consists of those involving resonance exchanges between photon pairs [131, 132, 133, 134]. 
Although here also the results depend on the models used, there is a constant feature that emerges from 
all the analyses that have been done: the contribution coming from the exchange of the pseudoscalars, 
7T°, r\ and r\' gives practically the final result. Other contributions [charged pion loops, vector, scalar, 
and axial resonances,...] are not only smaller, but also tend to cancel among themselves. 




Figure 11: Some individual contributions to hadronic light-by- light scattering: the neutral pion ex- 
change and the charged pion loop. There are other contributions, not shown here. 



Some of the results obtained for ' have been collected in Table 4. Leaving aside the first result 
[130, 3] shown there, which is affected by a bad numerical convergence [131], one notices that the sign 
of this contribution has changed twice. The first change resulted from a mistake in Ref. [131], that 
was corrected for in [133]. The minus sign that resulted was confirmed by an independent calculation, 
using the ENJL model, in Ref. [132]. A subsequent reanalysis [134] gave additional support to a 
negative result, while also getting better agreement with the value of Ref. [132]. 

Needless to say, these evaluations are based on heavy numerical work, which has the drawback of 
making the final results rather opaque to an intuitive understanding of the physics behind them. 
We [135] therefore decided to improve things on the analytical side, in order to achieve a better 
understanding of the relevant features that led to the previous results. Taking advantage of the 
observation that the pion exchange contribution a!jt' LxL,7r ' was found to dominate the final values 
obtained for af?' LxL ^ we concentrated our efforts on that part, that I shall now describe in greater 
detail. For a detailed account on how the other contributions to ajj 1 ' LxL ) arise, I refer the reader to 
the original works [131, 132, 133, 134]. 
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Table 4: Various evaluations of aft' LxL ' > x 10 11 and of the pion pole contribution aft' 



LxL;7i-°) 



x 10 11 . 



-260(100) 


constituent quark loop 


[130, 3] 


+60(4) 


constituent quark loop 


" [131] 


+49(5) 


7r ± loop, 7T° and resonance exchanges, aft' LxL,7r ) = +65(6) 


[131] 


-92(32) 


ENJL, aft' LxL ^°+^') = _85(13) 


[132] 


-52(18) 


tt^ loop, 7T° and resonance exchanges, and quark loop aft' LxL,7r ) = —55.60(3) 


[133] 


-79.2(15.4) 


ir^ loop, 7T° pole and quark loop, aft' LxL,7r ) = —55.60(3) 


[134] 


+83(12) 


7T°, rj and r/' exchanges only 


[135] 


+89.6(15.4) 


tt^ loop, 7T° pole and quark loop, aft' LxL,7r ** = +55.60(3) 


[136] 


+83(32) 


ENJL, aft' LxL ' 7r °+^') = +85(13) 


[137] 






Figure 12: The pion-pole contributions to light-by-light scattering. The shaded blobs represent the 
form factor J r 7r o 7 * 7 * . The first and second graphs give rise to identical contributions, involving the func- 
tion Ti(qi, q 2 ;p) in Eq. (5.30), whereas the third graph gives the contribution involving T 2 (qi,q2',p)- 



TT<y') i \ 



The contributions to 11^^(91,92,93) due to single neutral pion exchanges, see Fig. 12, read 

Zlivcxfi 9i <?2 e \ P ot 93 (91 + Q2) T 
£ppa(3 (92 + 93 ) /3 £v\*t 92<?3 
£p\af3 9l Q3 £ vp°T 92 (?1 + IsV ■ 



°7*7* {Ql + 92 + 93) 2 ) 

i'"^" 1 -" 2 -"*' " 1 (91 + 92) 2 - Ml 

,. ^V>7' 7 * (qI, (Q1 +12+ 93) 2 ) ^'7*7* («£> 9s) 



+^ 
+ i 



(92 + 9s) 2 - M2 
. ^07*7* {qi,qj) T n o rr (qj, (91 + 92 + 93) 2 ) 



(91 + 93) 2 - M2 



(5.28) 
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The form factor ./vVt*^!' ^2)) which corresponds to the shaded blobs in Fig. 12, is defined as 

i J d 4 xe^ x (n\T{j^x)j u (0)}\n°(p)) = e^fp? Fj> rr (q 2 , (p - q) 2 ) , (5.29) 

with ^ r 7r o 7 * 7 *(gi,<?2) = •^ r 7r°7*7* (^l) Inserting the expression (5.28) into (5.26) and computing the 
corresponding Dirac traces in Eq. (5.27), we obtain 



Q (h. LxL ;7 r°) = g 6 



d 4 <?i /" d 4 g 2 1 



7 (2vr) 4 7 (2vr) 4 g ^(9i+«2) 2 [(p + gi) 2 -m 2 ][(p-g 2 ) 2 -m 2 ] 
- / 1 ('/i •'/!•:/') 



•?V>7*7* (g l>(gl +g2) 2 ) ^ VU 7 » 7 » (?2) °) 

<Z 2 -M 2 

. 7 *(<? 2 ,gl) Jv> 7 . 7 .((gi + g 2 ) 2 ,0) 

H 7 : To Tin. 1 2{qi,q2;p) 

[qi + q2) 2 ~ 



(5.30) 



where T\ {q\ ,q2',p) and T 2 (q± ,q2',p) denote two polynomials in the invariants p ■ q\ , p ■ q2 , qi • q2 ■ 
Their expressions can be found in Ref. [135]. The former arises from the two first diagrams shown 
in Fig. 12, which give identical contributions, while the latter corresponds to the third diagram on 
this same figure. At this stage, it should also be pointed out that the expression (5.28) does not, 
strictly speaking, represent the contribution arising from the pion pole only. The latter would require 
that the numerators in (5.28) be evaluated at the values of the momenta that correspond to the pole 
indicated by the corresponding denominators. For instance, the numerator of the term proportional 
to T 1 (q 1 ,q 2 ;p) in Eq. (5.30) should rather read J^o^, (<??, (q\ + 2q x ■ q 2 + M 2 ) (M 2 , 0) with 

q\ = M 2 . However, Eq. (5.30) corresponds to what previous authors have called the pion pole 
contribution, and which they had found to dominate the final result. 

From here on, information on the form factor T^o^j* (q 2 , q 2 ) is required in order to proceed. The 
simplest model for the form factor follows from the Wess-Zumino-Witten (WZW) term [138, 139] that 
describes the Adler-Bell-Jackiw anomaly [140, 141] in chiral perturbation theory. Since in this case 
the form factor is constant, one needs an ultraviolet cutoff, at least in the contribution to Eq. (5.30) 
involving Ti, the one involving T 2 gives a finite result even for a constant form factor [131]. Therefore, 
this model cannot be used for a reliable estimate, but at best serves only illustrative purposes in 
the present context. 22 Previous calculations [131, 133, 134] have also used the usual vector meson 
dominance form factor [see also Ref. [143]]. The expressions for the form factor JF 7r o 7 » 7 * based on the 
ENJL model that have been used in Ref. [132] do not allow a straightforward analytical calculation 
of the loop integrals. However, compared with the results obtained in Refs. [131, 133, 134], the 
corresponding numerical estimates are rather close to the VMD case [within the error attributed to 
the model dependence]. Finally, representations of the form factor J 7 ^^*^*, based on the large- Nc 
approximation to QCD and that takes into account constraints from chiral symmetry at low energies, 
and from the operator product expansion at short distances, have been discussed in Ref. [144] . 
They involve either one vector resonance [lowest meson dominance, LMD] or two vector resonances 
[LMD+V], see [144] for details. The four types of form factors just mentioned can be written in the 

22 In the context of an effective field theory approach, the pion pole with WZW vertices represents a chirally suppressed, 
but large- Nc dominant contribution, whereas the charged pion loop is dominant in the chiral expansion, but suppressed 
in the large- N c limit [142] . 
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form [F n is the pion decay constant] 



T / 2 2 \ 7v 

J~ 7r o 1 * 1 *{q 1 ,q 2 ) — — 



^2 - M?, 



(5.31) 



For the VMD and LMD form factors, the sum in Eq. (5.31) reduces to a single term, and the corre- 
sponding function is denoted gM v (q 2 )- It depends on the mass My of the vector resonance, which will 
be identified with the mass of the p meson. For our present purposes, it is enough to consider only 
these two last cases, along with the constant WZW form factor. The corresponding functions f{q 2 ) 
and gM v (q 2 ) are displayed in Table 5. 

Table 5: The functions f(q 2 ) and gM v (q 2 ) of Eq. (5.31) for the different form factors. Nc is the 
number of colours, taken equal to 3, and F n = 92.4 MeV is the pion decay constant. Furthermore, 

_ N C K 



f(q 2 



9M V (q 



WZW 

VMD 

LMD 



N c 
4vr 2 F2 



Ml 



M v 



N c 

4vr 2 F 2 g 2 _ M 2 
q 2 + M y - C V 

q 2 ~ Ml 



We may now come back to Eq. (5.30). With a representation of the form (5.31), the angular integra- 
tions can be performed, using for instance standard Gegenbauer polynomial techniques [hyperspherical 
approach], see Refs. [145, 146, 76]. This leads to a two dimensional integral representation: 



a 



(h. LxL;7T°) 



,(^ ;i) - 



>°;2) 



0r°;i) + a (^;2) 



~Ll 



oo roo 

dQi / dQ 2 
o Jo 



roc roc 

= / dQi dQ 2 
Jo Jo 



w h (Qi,Q 2 ) f (1) (QlQl) 
+ w gi (M v ,Q 1 ,Q 2 ) g%{QlQ 2 2 ) 

w g2 (M,Q 1 ,Q 2 ) g { S{QlQl) 

M=M w ,Mv 



(5.32) 



(5.33) 
(5.34) 



The functions /^(Q?, Q 2 ), g$ v {Qj, Q 2 ), g^iQj, Q 2 ) and g$ {Q\ , Q 2 ) are expressed in terms of the 
functions given in Table 5, see Ref. [135], where the universal [for the class of form factors that have a 
representation of the type shown in Eq. (5.31)] weight functions w^, w gi , and w g2 in Eqs. (5.33) and 
(5.34) can also be found. The latter are shown in Fig. 13. 
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1 r : 0.06 r 



0.04- 




Figure 13: The weight functions appearing in Eqs. (5.33) and (5.34). Note the different ranges 
of Qi in the subplots. The functions and w gi are positive definite and peaked in the region 
Qi ~ Q2 ~ 0.5 GeV. Note, however, the tail in Wf 1 in the Qi-direction for Q2 ~ 0.2 GeV. The 
functions w g2 (M 7T ,Qi,Q2) and w g2 (Mv,Qi,Q2) take both signs, but their magnitudes remain small 
as compared to Wf 1 (Qx, Q2) and w gi (Mv, Qi, Q2)- We have used My = M p = 770 MeV. 



The functions Wf 1 and w gi are positive and concentrated around momenta of the order of 0.5 GeV. 
This feature was already observed numerically in Ref. [132] by varying the upper bound of the integrals 
[an analogous analysis is contained in Ref. [133]]. Note, however, the tail in Wfc in the Q\ direction 
for Q2 ~ 0.2 GeV. On the other hand, the function w 92 has positive and negative contributions in 
that region, which will lead to a strong cancellation in the corresponding integrals, provided they are 
multiplied by a positive function composed of the form factors [see the numerical results below]. As 
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can be seen from the plots, and checked analytically, the weight functions vanish for small momenta. 
Therefore, the integrals are infrared finite. The behaviours of the weight functions for large values of 
Qi and/or Q2 can also be worked out analytically. From these, one can deduce that in the case of 

the WZW form factor, the corresponding, divergent, integral for aft behaves, as a function of the 

ultraviolet cut off A, as aft°' 1] ~ C In 2 A, with [135] 

c = 3 (^) 2 (f) 2 = a0248 - (5 - 35) 

The log-squared behaviour follows from the general structure of the integral (5.33) for aft in the 
case of a constant form factor, as pointed out in [6]. The expression (5.35) of the coefficient C has 
been derived independently, in Ref. [147], through a renormalization group argument in the effective 
theory framework. 



[7T 1) 7T 2) 

Table 6: Results for the terms , and for the pion exchange contribution to the anomalous 

magnetic moment aft' LxL ' w ) according to Eq. (5.32) for the different form factors considered. In 
the WZW model, a cutoff of 1 GeV was used in the first contribution, whereas the second term is 
ultraviolet finite. 



Form factor 


(t u ;i) 


(t u ;2) 
ali 


aft' LxL;7rU) x 10 10 


WZW 


0.095 


0.0020 


12.2 


VMD 


0.044 


0.0013 


5.6 


LMD 


0.057 


0.0014 


7.3 



In the case of the other form factors, the integration over Qi and Q2 is finite and can now be performed 
numerically. 23 Furthermore, since both the VMD and LMD model tend to the WZW constant form 
factor as My — > 00, the results for aft in these models should scale as C In 2 My for a large resonance 
mass. This has been checked numerically, and the value of the coefficient C obtained that way is in 
perfect agreement with the value given in Eq. (5.35). The results of the integration over Q\ and Q2 
are displayed in Table 6. They definitely show a sign difference when compared to those obtained in 
Refs. [131, 133, 134, 143], although in absolute value the numbers agree perfectly. After the results 
of Table 6 were made public [135], previous authors checked their calculations and, after some time, 
discovered that they had made a sign mistake at some stage [136, 137]. The results presented in Table 
6 and in Refs. [135, 147] have also received independent confirmations [149, 148]. 

The analysis of [135] leads to the following estimates 

aft- LxL; " 0) = 5.8(1.0) x 10~ 10 , (5.36) 

and 

fl (h. LxL;7r°) = 5J x 1Q -14 (5>37) 
23 In the case of the VMD form factor, an analytical result is now also available [148]. 
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Taking into account the other contributions computed by previous authors, and adopting a conserva- 
tive attitude towards the error to be ascribed to their model dependences, the total contribution to 
coming from the hadronic light-by-light scattering diagrams amounts to 

af LxL) = 8(4) x 1(T 10 . (5.38) 



As a last remark, let me point out that the contribution depicted in Fig. 9 also involves the four point 
function H^ u \p{<li, Q2,Q3)- It would be interesting to have an evaluation of this contribution to ajt' Y ' p ' ^ 
based on the same models that were used to evaluate a^' LxL . The corresponding contribution arising 
from the neutral pion exchange has been evaluated in Ref. [148], but it is not obvious that it also 
dominates the complete result, since the kinematical configuration is different. This evaluation would 
also allow a direct comparison with the evaluations of Fig. 9 based on data. 



5.3 Electroweak contributions to a M 

Electroweak corrections to have been considered at the one and two loop levels. The one loop 
contributions, shown in Fig. 14, have been worked out some time ago, and read [150]-[154] 




Figure 14: One loop weak interaction contributions to the anomalous magnetic moment. 



Numerically, with G F = 1.16639(1) x 10~ 5 GeV~ 2 and sin 2 9 W = 0.224, 

a w(l) = 1Q4 8 x 1Q -ll ; (5>40) 

It is convenient to separate the two-loop electroweak contributions into two sets of Feynman graphs: 
those which contain closed fermion loops, which are denoted by a^ W ^' f , and the others, a^ W ^' b . In 
this notation, the electroweak contribution to the muon anomalous magnetic moment is 



a f = Q W(1) + a EW(2);f + a EW(2);b ^ 

I shall review the calculation of the two-loop contributions separately. 
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5.3.1 Two loop bosonic contributions 

The leading logarithmic terms of the two-loop electroweak bosonic corrections have been extracted 
using asymptotic expansion techniques, see e.g. Ref. [155]. In the approximation where sin 2 9\y — ► 
and Mh ~ M\y these calculations simplify considerably and one obtains 



EW(2);b _ "V « v 

" " V2 8vr 2 7T 



m 



(5.42) 



In fact, these contributions have now been evaluated analytically, in a systematic expansion in powers 
of sin 2 6w, up to C[(sin 2 9w) 3 ] > where In terms, In -^f- terms, ^ffi- In ^f- terms, terms and 

"V W H W H 

constant terms are kept [97]. Using sin 2 6 W = 0.224 and M H = 250 GeV , the authors of Ref. [97] find 



EW(2);b _ « v 

" ~ V2 8tt 2 vr 

G F ml 
y/2 8tt 2 \tt 



-5.96 In 



M, 



w 



mf 



+ 0.19 



n 



x [-78.9] = -21.1 x 10~ n , 
showing, in retrospect, that the simple approximation in Eq. (5.42) is rather good. 



(5.43) 



5.3.2 Two loop fermionic contributions 

The discussion of the two-loop electroweak fermionic corrections is more delicate. First, it contains 
a hadronic contribution. Next, because of the cancellation between lepton loops and quark loops in 
the electroweak £7(1) anomaly, one cannot separate hadronic effects from leptonic effects any longer. 
In fact, as discussed in Refs. [156, 157], it is this cancellation which eliminates some of the large 
logarithms which, incorrectly were kept in Ref. [158]. It is therefore appropriate to separate the 
two-loop electroweak fermionic corrections into two classes: One is the class arising from Feynman 
diagrams containing a lepton or a quark loop, with the external photon, a virtual photon and a virtual 
Z° attached to it, see Fig. 15. 24 The quark loop of course again represents non perturbative hadronic 
contributions which have to be evaluated using some model. This first class is denoted by a^ W ^'\i; q). 
It involves the QCD correlation function 

W^ p (q,k) = J d A xe^J d 4 ye^-^y(n\T{j,(x)Ai z \y)j p (0)}\n), (5.44) 

with k the incoming external photon four-momentum associated with the classical external magnetic 

(z) 

field. As previously, j p denotes the hadronic part of the electromagnetic current, and A p is the axial 
component of the current which couples the quarks to the Z° gauge boson. The other class is defined 
by the rest of the diagrams, where quark loops and lepton loops can be treated separately, and is 
called a^ W ^' f (residual). 



If one works in a renormalizable gauge, the contributions where the Z is replaced by the neutral unphysical Higgs 
should also be included. The final result does not depend on the gauge fixing parameter £z, if one works in the class of 
't Hooft gauges. 
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function W^ up (q, k). 



The contribution from a^ W ^' f (residual) brings in factors of the ratio rn^/My^. It has been evaluated, 
to a very good approximation, in Ref. [157], with the result 



1 



2 sin 2 9w 



5 m? , m 



+ A 



Higgs 



(5.45) 



where A Higgs denotes the contribution from diagrams with Higgs lines, which the authors of Ref. [157] 
estimate to be 

A Higgs = -5.5 ±3.7, (5.46) 

and therefore, 



a EW(2);f (residual) = x [-21(4)] = -5.6(1.4) x KT 11 



(5.47) 



Let us finally discuss the contributions to a^ W ^' f (.£; q). Here, it is convenient to treat the contributions 
from the three generations separately. The contribution from the third generation can be calculated 
in a straightforward way, with the result [156, 157] 



a 
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-ii 



M| ln ml 

w? t Ml i 



(5.48) 



In fact the terms of O 



^4 

mj Mi 



and O 



have also been calculated in Ref. [157]. There are 

in principle QCD perturbative corrections to this estimate, which have not been calculated, but the 
result in Eq. (5.48) is good enough for the accuracy required at present. The contributions of the 
remaining charged standard model fermions involve the light quarks u and d, as well as the second 
generation s quark, for which non perturbative effects tied to the spontaneous breaking of chiral 
symmetry are important [156, 159]. The contributions from the first and second generation are thus 
most conveniently taken together, with the result 
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, Ml A , Ml 11 8 2 
-31n^+41n— f - — + -vr 2 
m^ m~ 6 9 

m 2 in M| N 



M| ml 



+4.67(1.80) +0.04(2) I 
-28.5(1.8)] = -7.6(5) x 10~ n , (5.49) 



where the first line shows the result from the e loop and the second line the result from the ll loop 
and the c quark, which is treated as a heavy quark. The term between brackets in the third line is 
the one induced by the anomalous term in the hadronic three point function W^ up (q,k) The other 
contributions have been estimated on the basis of an approximation to the large-Afc limit of QCD, 
similar to the one discussed for the two-point function H(Q 2 ) after Eq. (5.15), see Ref. [159] for 
details. 

The result in Eq. (5.49) for the contribution from the first and second generations of quarks and 
leptons is conceptually rather different from the corresponding one proposed in Ref. [157] , 
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EW(2) ;ff d x = GV r?V a 



o, M z „, M z , M z 5 
-3 In —| + 41n^ - In— f - - 



Ml Ml Ml 11 8 , 

-31n^ + 41n^-ln^- — + -vr 2 -6 
m^ m~ mi. b 9 

= % 2L x [-31.9] = -8.5 x 10- 11 , (5.50) 

y/2 OTT Z IT 

where the light quarks were, arbitrarily, treated the same way as heavy quarks, with m u = md = 
0.3 GeV , and m s = 0.5 GeV . Numerically, the two expression lead to similar results, though. A more 
recent analysis [160] provides a non perturbative treatment of the light quark sector, and gives 

a™^(e, ^ u, d, s, c) = ^ ^ H x [-24.6] = -6.6 x lO" 11 . (5.51) 

The difference between the two results in Eqs. (5.49) and (5.51), which is numerically very small, is 
connected to interesting issues involving the anomalous (VVA) three point function in QCD. 25 

The authors of Ref. [160] have also performed a detailed renormalization group analysis of the leading 
logarithm contributions at three loops 26 , and found them to be negligible. Taking into account other 
small effects that were previously neglected, their final value reads [160] 

a|f w = 15.4(3) x 10- 10 , (5.52) 

which shows that the two-loop correction represents indeed a reduction of the one-loop result by an 
amount of 23%. The final error here includes uncertainties in the hadronic part, the variation of the 
Higgs mass in a range 114 GeV < Mu < 250 GeV, the uncertainty on the mass of the top quark, and 
unknown three loop effects. 

25 Actually, the discussion centers around the transverse, i.e. non anomalous part of this QCD correlator, and is related 
to the existence of non renormalization theorems for it [161, 160, 162]. 
26 See also [163]. 



36 



5.4 Comparison with experiment 



We may now put all the pieces together and obtain the value for a„ predicted by the standard model. 
We have seen that in the case of the hadronic vacuum polarization contributions, the latest evaluation 
[122] shows a discrepancy between the value obtained exclusively from e + e~ data and the value that 
arises if r data are also included. This gives us the two possibilities 

a^ M (e + e~) = (11 659 167.5 ± 7.5 ± 4.0 ± 0.4) x 1(T 10 
oJ m (t) = (11 659 192.7 ± 5.9 ± 4.0 ± 0.4) x 10^ 10 , 

(5.53) 

where the first error comes from hadronic vacuum polarization, the second from hadronic light-by- 
light scattering, and the last from the QED and weak corrections. When compared to the present 
experimental average 

a e ^> = (11 659 203 ±8) x 10~ 10 (5.54) 

there results a difference, 

< P " a^ M ( e+e ~) = 35.5(11.7) x KT 10 , 
a cx P _ a SM( r ) = 10 .3(io.7) x HP 10 , 

which represents 3.0 and 1.0 standard deviations, respectively. 

Although experiment and theory have now both reached the same level of accuracy, ~ ±8 x 10~ 10 or 
0.7 ppm, the present discrepancy between the e + e~ and r based evaluations makes the interpretation 
of the above results a delicate issue as far as evidence for new physics is concerned. Other evaluations 
of comparable accuracy [119, 120, 58] cover a similar range of variation in the difference between 
experiment and theory. Furthermore, the value obtained for a^ M (e + e~) relies strongly on the low- 
energy data obtained by the CMD-2 experiment, with none of the older data able to check them at 
the same level of precision. There seems to be an error in these data from the CMD-2 experiment [8], 
but this clearly needs to be confirmed. In this respect, the prospects for additional high statistics data 
in the future, either from KLOE or from BaBar, are most welcome. On the other hand, if the present 
discrepancy in the evaluations of the hadronic vacuum polarization finds a solution in the future, and 
if the experimental error is further reduced, by, say, a factor of two, then the theoretical uncertainty 
on the hadronic light-by-light scattering will constitute the next serious limitation on the theoretical 
side. It is certainly worthwhile to devote further efforts to a better understanding of this contribution, 
for instance by finding ways to feed more constraints with a direct link to QCD into the descriptions 
of the four-point function U^ upcT (qi,q2,q3)- 



6 Concluding remarks 

With this review, I hope to have convinced the reader that the subject of the anomalous magnetic 
moments of the electron and of the muon is an exciting and fascinating topic. It provides a good 
example of mutual stimulation and strong interplay between experiment and theory. 
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The anomalous magnetic moment of the electron constitutes a very stringent test of QED and of 
the practical working of the framework of perturbatively renormalized quantum field theory at higher 
orders. It tests the validity of QED at very short distances, and provides at present the best determi- 
nation of the fine structure constant. 

The anomalous magnetic moment of the muon represents the best compromise between sensitivity to 
new degrees of freedom describing physics beyond the standard model and experimental feasibility. 
Important progress has been achieved on the experimental side during the last couple of years, with 
the results of the E821 collaboration at BNL. The experimental value of is now known with an 
accuracy of 0.7ppm. Hopefully, the Brookhaven experiment will be given the opportunity to reach its 
initial goal of achieving a measurement at the 0.35 ppm level. 

As can be inferred from the examples mentioned in this text, the subject constitutes, from a theoretical 
point of view, a difficult and error prone topic, due to the technical difficulties encountered in the 
higher loop calculations. The theoretical predictions have reached a precision comparable to the 
experimental one, but unfortunately there appears a discrepancy between the most recent evaluations 
of the hadronic vacuum polarization according to whether r data are considered or not. Hopefully, this 
situation will be clarified soon. Hadronic contributions, especially from vacuum polarization and from 
light-by-light scattering, are responsible for the bulk part of the final uncertainty in the theoretical 
value a^ M . Further efforts are needed in order to bring these aspects under better control. 
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